Galaxy Formation: Warm Dark Matter, Missing Satellites, 
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Abstract. We present warm dark matter (WDM) as a possible solution to the 
missing satellites and angular momentum problem in galaxy formation and introduce 
improved initial conditions for numerical simulations of WDM models, which avoid 
the formation of unphysical haloes found in earlier simulations. There is a hint, that 
, the mass function of satellite haloes has been overestimated so far, pointing to higher 

values for the WDM particle mass. 

<: 

The cold dark matter (CDM) scenario for structure formation is well 
CN ■ established, and is particularly successful in explaining the formation 

, of large-scale structure and galaxies. But there exist problems on small 

' (i.e. galactic) scales, among them: (1) The missing satellites problem: 

■ CDM produces too many small galaxies, e.g. about five times as many 
satellites as are observed in the local group (Klypin et al., 1999; Moore 

■ et al., 1999; Kamionkowski and Liddle, 2000). (2) The angular momen- 
. turn problem: Galaxies in CDM simulations consistently have smaller 
\ specific angular momenta, and therefore have smaller disks than what 

O is observed (Sommer-Larsen and Dolgov, 2001 and references therein). 

q \ One possible solution to these problems is to go from CDM to warm 

dark matter (WDM) . The free-streaming motion of the WDM particle 
reduces power on small scales, but keeps it unchanged at long wave- 
lengths, not disturbing the predictions of CDM for the formation of 
large-scale structure. This leads to the formation of fewer low-mass 
^ , systems, explaining the missing satellites problem, and to fewer merg- 

5h ' ing events, during which subclumps, which are assembled into galactic 

disks, would loose orbital angular momentum and energy by dynam- 
ical friction. Typically, WDM particle masses of the order l/i 5//4 keV 
(corresponding to a free-streaming mass of Mf ~ 3 • lO 11 ft, _1 M0 for 
S7o = 0.3) are necessary to explain the missing satellites (Kamionkowski 
and Liddle, 2000; Colin et al., 2000; White and Croft, 2000; Bode et 
al., 2001; Governato et al., 2002) and to reduce the discrepancy of the 
specific angular momenta from a factor of 10 or more to about a factor 
of 2 (Sommer-Larsen and Dolgov, 2001). 

Thus to check WDM against observations, it is necessary to study, 
how many small haloes there are, and how the large disk galaxies 
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Figure 1. Dark matter particles at redshifts z = and z = 2 within a slice of size 
10 x 10 x 0.325 (/i _1 Mpc) for two AWDM simulations starting with identical phases 
from grid and glass initial conditions, respectively. The regularly spaced, unphysical 
haloes along filaments in grid simulations are clearly visible at z — (in particular 
along the filament running from the halo in the center to the one close to the lower 
right corner), and how they form as artifacts due to trains of particles falling onto 
the forming filaments at high redshifts. They are not present in the run starting 
from a glass. 

are formed. This requires numerical simulations, even for the rela- 
tively simple task of calculating the halo mass function in a WDM 
scenario. There, structure formation does not follow the hierarchical 
picture well-known from CDM, since haloes smaller than the free- 
streaming scale form later than the larger ones by non-linear transfer of 
power from large to small scales. That means that analytical schemes 
(e.g. Press-Schechter theory and its variations) can not be applied to 
WDM. Figure 1 in Gotz and Sommer-Larsen (2002) clearly shows the 
discrepancy between Press-Schechter theory and the numerical mass 
function. 

But numerical simulations with WDM have problems of their own. 
To check the effect of reduced power on small scales, they have to be 
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Figure 2. The cumulative number density of subhaloes as a function of their maxi- 
mum circular velocity V^ max . The subhaloes have to be located within AQ0h~ kpc of 
the density peak of the host halo, which are selected to be Milky Way-like with K.max 
between 220 and 280 km/s. The solid and dotted lines show the numerical results 
for the ACDM and AWDM (with a free streaming mass of M f = 1.5 ■ lO^^Me) 
simulations with glass initial conditions presented in this paper, together with their 
1<7 Poisson errors. The gray area indicates, where these simulations are incomplete. 
The observed number density for the Milky Way and M31 from Klypin et al. (1999) 
is shown by filled circles, and their numerical result for a ACDM-model by the 
dashed line, which agrees with our result. 

set up such, that the mean particle separation is (much) smaller than 
the wavelength, below which the power spectrum is cut off. Thus initial 
displacements and velocities of neighboring particles are highly corre- 
lated. Since at high redshifts, matter moves perpendicularly onto the 
forming filaments, there will be "trains" of particles falling onto them 
in the usual set-up, where one starts from a regular cubic mesh. They 
create unphysical haloes along the filaments, with a separation given by 
the grid spacing projected onto the filament (Gotz and Sommer-Larsen, 
2002). These "beads on a string", and how they form, can clearly be 
seen in the left hand panels of Figure 1. To avoid this problem, the grid 
structure needs to be broken by starting from glass-like initial condi- 
tions (White, 1996), where the particles are irregularly distributed, but 
still (almost) evenly spaced. The right hand panels in Figure 1 show 
the now much smoother distribution of matter along the filaments. 

To calculate a more reliable satellite mass function, not affected by 
the presence of unphysical haloes, we ran a WDM simulation start- 
ing with glass initial conditions. We chose a A-cosmology (VLq = 0.3, 
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J7a = 0.7, h = 0.65) with a cluster-abundance normalized power spec- 
trum (o"8 = 1.0) and a free streaming mass of Mf = 1.5 • lO 11 /?^ 1 M@, 
corresponding to a WDM particle mass of 1.2/i 5 / 4 keV close to the value 
found by other authors. We used the publicly available Hydra code 
(Couchman et al., 1995) with 128 3 dark matter particles in a 10/t _1 Mpc 
box. A corresponding CDM simulation was also run. Candidate host 
haloes were identified with friends-of- friends, and Milky Way-like hosts 
with a maximum circular velocity K,max between 220 and 280 km/s 
were selected — 8 in the AWDM and 9 in the ACDM run. Within these 
haloes, subhaloes were found with a variation of the HOP algorithm 
(Eisenstein and Hut, 1998), and Figure 2 shows the cumulative number 
densities of those satellites, which lie within 400/i -1 kpc of the density 
peak of the host halo, cts ci function of their max* 

As expected, CDM 

seems to overpredict the number of satellites compared to the observed 
distribution for the Milky Way and M31 (Klypin et al., 1999), but our 
improved WDM simulations without the unphysical haloes now give 
numbers which are too low. This hints a lower Mf and higher WDM 
particle mass than has been suggested so far. But that result has to be 
taken with caution, since the identification of subhaloes is non-trivial, 
and we have not yet checked the effects of different algorithms on it. 

In conclusion, WDM can solve the missing satellites and angular 
momentum problems, but improved initial conditions starting with a 
glass-like distribution of the particles are necessary for numerical stud- 
ies to avoid the appearance of unphysical haloes which are purely grid 
artifacts. Thus there is the possibility that the satellite mass function 
has been overestimated in previous studies, and that the WDM particle 
mass actually has to be higher. But it should be mentioned that there 
are other open questions for WDM. E.g., high velocity clouds could 
be the missing satellites (Braun and Burton, 1999), CDM with feed- 
back from star formation may solve the angular momentum problem 
(Sommer-Larsen et al., 2002; Thacker and Couchman, 2001), and the 
late structure formation in WDM can be a problem, if the trend of 
finding quasars at ever higher redshifts continues. 
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